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Abstract. We revisit the problem of a two-dimensional polymer ring subject to an inflating pressure dif- 
ferential. The ring is modeled as a freely jointed closed chain of TV monomers. Using a Flory argument, 
mean-field calculation and Monte Carlo simulations, we show that at a critical pressure, p c ~ TV -1 , the 
ring undergoes a second-order phase transition from a crumpled, random-walk state, where its mean area 
scales as (A) ~ TV, to a smooth state with (A) ~ TV 2 . The transition belongs to the mean-field universality 
class. At the critical point a new state of polymer statistics is found, in which (A) ~ TV 3//2 . For p > p c we 
use a transfer-matrix calculation to derive exact expressions for the properties of the smooth state. 

PACS. 36.20.Ey Macromolecules and polymer molecules: Conformation (statistics and dynamics) - 
05.40.Fb Random walks and Levy flights - 64.60.-i General studies of phase transitions 



1 Introduction 

Considerable theoretical efforts were directed during the 
1980s and 1990s at random polymer rings constrained to 
a plane, both as a fundamental problem of statistical me- 
chanics [lB2] |3 ]|4l | $j and as a highly idealized model for 
membrane vesicles |(3,7 8,9 1 0lllllll2llT 3 14 . In addition, 
there was much interest in transitions between crumpled 
and smooth states of membranes |15U16j . 

Two bodies of work, in particular, concerned pressur- 
ized two-dimensional (2D) rings. Fisher et al. studied both 
inflated and deflated, closed, 2D self-avoiding walks using 
Monte Carlo (MC) simulations and scaling analysis 013 
For a finite ring of TV monomers, subject to an inflating 
pressure differential p > 0, three regimes were found: (i) 
a weak-inflation crumpled regime for < pN 2v < 1 (y = 

■ 3/4 being the swelling exponent of 2D self-avoiding walks), 
where the mean area scales as (A) ~ N 2v Ja(pN 2u ); (ii) a 

, strong- inflation crumpled regime for 1 < pN 2v <C N 2v ~ l , 
where the same scaling holds but the scaling function 
/a becomes a power law; and (iii) a smooth regime for 
pN 2u > TV 2 " -1 , where (A) — TV 2 . The crossovers between 
these regimes are gradual with no phase transitions. In the 
thermodynamic limit, defined as TV — > oo and p — ► such 
that pN 2v is finite, the range of the scaling regime (ii) in 
the case of self-avoiding walks (y — 3/4) becomes infinite. 
Note, however, that in the random-walk case [y = 1/2) 
this regime disappears; we will show below that for ran- 
dom walks regime (ii) turns, in fact, into a second-order 
phase transition between regimes (i) and (iii). 

The second body of works consists of analytical studies 
of closed, pressurized, 2D random walks, first presented by 
Rudnick and Gaspari [9"lll0llll| . In this model the poly- 



mer comprises a set of Gaussian springs of fixed elas- 
tic constant and no bending rigidity. At a critical pres- 
sure p c ~ TV -1 the mean area was found to diverge (i.e., 
the ring inflates to infinity). This divergence, obviously, 
is made possible by the extensibility of the chain, i.e., 
the ability of the springs in this model to be infinitely 
stretched. 

In the current work we revisit the Rudnick-Gaspari 
model while imposing inextensibility of the chain, as is 
actually appropriate for real polymers or vesicles. This 
constraint changes the infinite inflation at p = p c into a 
second-order phase transition between a crumpled, random- 
walk state and a smooth one. The model is defined in Sec- 
tionEI We begin the investigation in Section|3with a sim- 
ple Flory argument which, nevertheless, yields the correct 
physics. In Section 0] we present a mean-field calculation 
which accurately captures the behavior below and at the 
transition. In Section[3]we derive exact expressions for the 
behavior at p 3> p c using a transfer-matrix technique. The 
behavior found in Sections 13151 is verified in Section by 
MC simulations. Finally, the results are summarized and 
discussed in Section 



2 Model 

The system under consideration is illustrated in Fig.Q] A 
closed, two-dimensional, freely jointed chain of TV monomers 
is subject to an inflating 2D pressure differential p > 
between its interior and exterior. The monomers are con- 
nected by rigid links of length I. We define I = 1 as the 
unit length and the thermal energy U-qT as the unit en- 
ergy. (Thus, p is measured in units of k^,T /I 2 .) The chain 
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is ideal, i.e., there are no interactions between monomers. 
(Effects related to self-avoidance will be briefly discussed 
in Sections |21 and 0) fn addition, no bending rigidity is 
taken into account, i.e., the chain is freely jointed. A con- 
figuration of the ring is defined by a set of 2D vectors 
{?j}j=o...N specifying the positions of the monomers. The 
condition that the chain be closed is expressed by ro = rjf. 




Fig. 1. Schematic illustration of the ring and its parameters. 



The probability of a specific configuration is 

N 

P ({ rj },p) ex e^K'i" J] %i - r ^il - !)> 

i=i 

where A is the area enclosed by the ring. As in previous 
works 2.3.9 fTUllll| . we take A as the algebraic area rather 
than the geometrical one, 

1 N 

A[{r i }] = -^(r,-_ 1 xr i )-z, (2) 

3=1 

where z is a unit vector perpendicular to the plane of the 
ring. This area may take both positive and negative val- 
ues. At zero pressure both signs are equally favorable and 
the mean area must vanish. At high pressures the prob- 
ability of configurations with negative A is exponentially 
small in p\A\, and whether one takes the algebraic or ge- 
ometrical area will become statistically insignificant. (We 
shall further discuss this assumption in Section 0) 

Using Eqs. Q and J2J), we write the partition function 
of the ring as 

r N 

Z(p,N) = / Hdvje^-^^Sdvj-r^l-l), (3) 

J 3=1 

with ro — tn- 



3 Flory Argument 

We begin the analysis with a simple Flory argument ^7] 
which captures most of the physics to be more rigorously 
treated in the following sections. The free energy of the 
ring (in units of ksT) is expressed as a function of R, the 
radius of the statistical cloud of monomers (i.e., the root- 
mean-square radius of gyration). We divide it into three 
terms, 

F(R) = F e \ + Fjncxt + F p 
F el ~ R 2 /N 

-Finext ~ R /N 3 

F p ~ - P R 2 . (4) 

The elastic term, F c \, is the usual entropic-spring free en- 
ergy of a Gaussian chain |TZ]. The second term is the lead- 
ing non-Gaussian correction due to the inextensibility of 
the chain [see Appendix, Eq. (|A. 51) ] . needed here to sta- 
bilize the ring against infinite expansion. The last term is 
the pressure contribution, where the mean area of the ring 
is taken as proportional to R 2 

Equation J3J has the form of a Landau free energy, 
describing a second-order transition at p = p c ~ A -1 . 
Since the critical pressure depends so strongly on system 
size, we use hereafter the rescaled pressure p = p/p c ~ pN, 
and define the thermodynamic limit as N — > oo and p — > 
such that p is finite. For p < 1 -Fmext is negligible, and R 
has a Gaussian distribution with (R 2 ) ~ N(l — p)^ 1 . For 
p > 1 we have R 2 ~ N 2 (p — 1). Thus, defining an order 
parameter M — R/N, we find in the thermodynamic limit 

M ~{°(P~iy, p = i/2 Itl (5) 

At the critical point itself F = i*i ncxt ~ R 4 /N 3 , and R 
has a non-Gaussian distribution with 

(i? 2 G5=l))~A 2l % ^ c = 3/4. (6) 

Note that the competition between F c \ and F p , lead- 
ing to the second-order transition, is unique to 2D. In 
addition, when an excluded-volume term, F cv ~ N 2 /R 2 , 
is added to the free energy, Eq. Ql, the transition is re- 
moved. This agrees with previous studies of self-avoiding 
rings [UIZIIHjj which did not report any phase transition 
upon increasing pressure. 



4 Mean-Field Theory 

In this section we calculate the partition function of the 
freely jointed ring, Eq. 10, by relaxing the rigid delta- 
function constraints on link lengths into harmonic poten- 
tials. The spring constant A of the links is chosen so as 
to make the root-mean-square length of a link equal to 
1 = 1. This type of approximation, first suggested by Har- 
ris and Hearst [TH| , was successfully employed in studies of 
the Karatky-Porod worm- like-chain model |18U19j , where 
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it was shown to be equivalent to a mean-field assump- 
tion (for the field conjugate to the rigid link-length con- 
straints). The partition function contains now only Gaus- 
sian terms and, therefore, can be calculated exactly, 



r N 

Z(p,N,X) = J Y[dr d e? 



p(r.,_i xrj)-z— A(r 3 — r,_i) 

(7) 
1 Np 



A^4Asin(^)- 

(The spring constant A is in units of k^T/l 2 .) Details of 
the calculation can be found in Ref. ^Uj- This result can 
be obtained also by analogy to the quantum propagator 
of a charged particle in a magnetic field [23j- The mean 
area is obtained by differentiation with respect to p, 




(A(p,N,X)) = 



dlnZ 
dp 



1 

P 



N cot (*£) 



4A 



(8) 



For A = 1 Eq. JSJ is the same as the result obtained by 
Rudnick and Gaspari |9lll()llll| , exhibiting a divergence at 



Pc 



Att/N. 



(9) 



Yet, in our case A is not fixed but is to be determined 
self-consistently to ensure the softened inextensibility con- 
straint. It is clear that, as the pressure increases, the springs 
must become stiffer to satisfy this constraint. To impose X(p,N ^> 1) 
the constraint, we calculate the mean-square link length 
and set it to 1, 



Fig. 2. The mean area in units of A max ~ N 2 as a function of 
the rescaled pressure (in units of p c ~ iV _1 ) as obtained from 
the mean-field approximation. The dashed curves are calcu- 
lated numerically using Eqs. fllOP and Hill , whereas the solid 
curves (calculated only for p > 1) present the approximation 
given by Eq. 1151 . Calculations were performed for N — 10 
(thin curves) and N = 10 (thick curves). 



cot[ATp/(4A)] ~ [Np/(4X) - it}- 1 . This gives 



i 



(13) 



((r i -r J -_ 1 ) 2 ) = - 



1 d\nZ 




where p = p/p c — pN / (4-7r) is the rescaled pressure. In the 
thermodynamic limit Eq. I|13|l reduces to the continuous 
but nonanalytic function 



4A 



(10) 



1, 



X(p,N 



oo 




(14) 



thus obtaining a transcendental equation for X(p, N). Equa- 
tions (JHJ) and (|10|l are combined to yield a simpler expres- 
sion for (A) as a function of A, 



(A(p,N,X)) = 



N(X-l) 



(11) 



Equation l(DB*)l should be regarded as an asymptotic ex- 
pression in the limit Np/(AX) — > n, which turns out to 
be valid for any p <jt. p c - (A Taylor expansion around this 
point fails because of the nonanalyticity inferred above.) 
Substituting Eq. 1)13(1 in Eq. (|11J) yields an approximate 
expression for (A) as a function of p and N, 



Numerical solution of Eq. Hl()|l for A [in the range A > 
Np/(4ir)] and substitution of the result in Eq. (|11|) yield 
the mean area as a function of p and N. (See dashed curves 
in Fig. El) 

For very low pressures, p<p c , we expand Eq. JHJ to 
first order in p to get 



0-1)2 + £(0+1) 



1 

N- 



(15) 



4?r 



2f> 



(A(p^ Pc ,N)) = -N 2 p, 



(12) 



In the thermodynamic limit, the behavior of (A) around 
and above the critical pressure is obtained from Eq. (|15fl 



i.e., a linear dependence on p as expected from linear re- 
sponse. (Recall that at p — the mean algebraic area 
vanishes.) 

For higher pressure we obtain a good approximation 
for X(p, N) in the limit N » 1 by substituting in Eq. lfl"T)|l 



(A) = 



N 1 

47T p(l— p) 

in 

N 2 p-1 P_ 

4-7T p 



N 1 

47T 1— p 



1 -p > N- 1 ' 2 
\l-p\ « N- 1 ' 2 
p- 1> TV" 1 / 2 , 



(16) 
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revealing a continuous (second-order) transition. Below 
the transition we get the same behavior as in the Rudnick- 
Gaspari model PITH) . (A) ~ A(l - py 1 . Yet, due to 
the inextensibility in our model, the increase of (A) as 
p — ► l - breaks at |1 — p\ ~ TV -1 / 2 . In the transition 
region, |1 — p\ <C A -1 / 2 , which shrinks to a point in 
the thermodynamic limit, we find (A) ~ N 2Uc ,v c — 3/4. 
Above the transition the ring reaches a smooth state with 
(A) ~ A 2 (p — l)/p. All of these results agree with the 
findings of the Flory argument presented in Section[3]once 
we identify (A) ~ (Rg)- As p is increased to infinity, (A) 
tends, as it should, to 



(A(p — > oo, AO) = A n 



A 2 

47T ' 



(17) 



which is the area of a circle of perimeter A. 

Figure El shows the dependence of (A) on p for A = 10 
and 10 5 calculated both from the numerical solution of 
Eqs. i|10|) and (|llf) . and using the approximate expression 
H15(l . For large A the critical behavior becomes appar- 
ent, with a transition between two distinct states — one in 
which (A) ~ A and, hence, in units of A max ~ A 2 , the 
mean area vanishes for A — > oo, and another with a mean 
area proportional to A max ~ A 2 . As can be seen in Fig. 
12 the approximate expression (|15J) is practically indistin- 
guishable from the numerical solution for p > 1, even for 
small A. 

The compressibility (defined with respect to the re- 
duced pressure p) is obtained from Eq. i|TT)|) as 



1 OA 



. .... p-n=" 

Ifl* >b 



ir 7 , 7 = i- (is) 

At the critical point itself the compressibility diverges as 

(19) 



K{p=\)= l -N 1 ' 2 . 



To calculate the mean-square radius of gyration, (R 2 ) = 
A -1 r 2 (where Yj are measured with respect to the cen- 
ter of mass), we add a h r 2 term to the Hamiltonian of 
Eq. (J2J) and differentiate the resulting partition function 
with respect to h. This yields 



(R 2 g (p,N,X)) = 



1 d\nZ(p,N,\,h) 



A 



Oh 



h=0 



4A-Apcot(^) 
Np 2 



(20) 



For p< 1, combining this result with Eqs. I|12l) and (|14l) . 

we get (i? 2 ) = N/12, which is the well known result for the 
mean-square radius of gyration of a Gaussian ring (e.g., 
\17\). For large pressures, p > 1, we have from Eq. 11411 

A — — ^+ p — pN / (47r) , thereby recovering the relation for 
an average circle, 7r(i? 2 ) = (A). 

Figure shows the dependence of (R 2 ) on A at fixed 
pressure. The data were obtained by substituting the nu- 
merical solution for A [Eq. ljT3|) ] in Eq. (J2DJ. The scaling 
of (R 2 ) changes at the critical point A c = 47r/p. Below 



the critical point, A < A c , (R 2 ) ~ A, as in a Gaussian 
chain. Above it (R 2 ) ~ A 2 , as in a stretched chain. As 
p decreases, the transition becomes sharper. At exactly 
A = A c (i? 2 ) scales as A 3 / 2 . Thus, the analysis of (R 2 ) 
yields the same scaling with A as obtained for (A) . 
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Fig. 3. Mean-square radius of gyration as a function of JV at 
fixed values of p = 10 -5 , 10 -3 and 10 _1 (thick, dashed and thin 
curves, respectively) as obtained from the mean-field theory. 
For N < N c — 4%/p (marked with arrows) (R 2 S ) ~ N, while 
for N > N c (i?g) ~ N 2 . At the critical point (diamonds) (Rg) 
is proportional to N 3 ^ 2 . The dotted lines show, from bottom to 
top, the dependencies for an unpressurized ring ((Rg) = N/12), 
for a ring at the critical pressure [{Rg) = 7V 3//2 /(47r 2 )] and for 
a stretched circle [(R 2 S ) = N 2 /{4tt 2 )]° 



5 Transfer- Matrix Formulation 

Since the interactions between all adjacent monomers are 
identical, including the one between the first and the Ath, 
the partition function may be rewritten in the following 
transfer-matrix form: 



N 



Z(p, A) = Tr TT Tfo-i.rj), r = v N 

T(r,r') = e5 (rxr,) - 2 (5(|r-r'| - 1). 
Solution of the associated eigenvalue problem, 

l«<?(r) = J dr'T(r,r')<P(r'), 



(21) 



(22) 



(in particular, finding the two eigenvalues /i of largest ab- 
solute value) will yield the exact solution of the model. 

Two properties of the operator T are readily noticed: it 
is non-Hermitian, T(r, r') ^ T(r',r), and it is rotational- 
invariant. To exploit the invariance to rotations, we change 
to polar coordinates, r = (p, ip), and separate variables as 
&(r) — x(p) e * mV j where m = 0, ±1,±2, ... to maintain 
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periodicity in ip. Equation (|22|l can then be integrated 
over angles to give 



PX{p) 
*W) = 
A(p,p') = 



p+i 



dp'f(p,p') X (p') 



p-i 



2 cosh (ppp'^Jl- A 2 /2 



im cos 



A) 



pVl - A 2 



(23) 



P l + P 



2pp> 



thus reducing the original operator to the one-dimensional 
operator T. 

Unfortunately, we have not been able to diagonalize T 
for any p. For p 3> p c , nevertheless, the ring has stretched 
configurations, and we can assume that the distances of 
all monomers from the center of mass are much larger 
than the link length, /)> 1. We expand x(p') around p to 
zeroth order, change variables according to p' = p + sin 9, 
and integrate over 9 to get 

fi^x(p) = 2TrI (pp/2)x(p), (24) 

where Iq is the zeroth-order modified Bessel function of 
the first kind. Thus, at this order of approximation, the 

eigenfunctions have the form x^ (p) = [-^ r /(2 7r Pfc)] 1 ^ 2 ^(jO— 
Pk) with a continuous spectrum of eigenvalues. The spec- 
trum is bounded from above by /iSax = 2-n:Ia(pp max /2) = 
27rio[pAT/(47r)], where p max — N/(2n) is the radius of a 
perfect circle of perimeter N. (The value of /0 ma x can also 
be obtained from the condition that J \x(p)\ 2 dp > 1, i.e., 
that going from the center of mass outward, one must 
cross the ring at least once.) 

Within the zeroth-order approximation, and in the ther- 
modynamic limit, the partition function is given by 



Z^(p,N) = [p 



(0) iJV 
max J 



[2irI (pN/4ir)} 



N 



(25) 



This result has a straightforward interpretation. As shown 
in the Appendix [Eq. l|A.l(l ]. it is identical to the partition 
function of a 2D, open, freely jointed chain subject to a 
tensile force / = pN/(4ir). This force is just the tension 
associated with a Laplace pressure p acting on a circle of 
radius /9 max ~ N. The mean area is obtained from Eq. 
as 

d\nZ N 2 h{p) 



(A(p,N)) 



(26) 



dp 4tt I (p) ' 

which saturates, as expected, to A max — N 2 /(An) as p 
oo. The approach to saturation is given by 

1 

w 

Eq. (QiJ), by a 



(A(p » 1))/A D 



1 



(27) 



which corrects the mean-field prediction, 
factor of 2. 



of gyration (R 2 ) and mean-square area fluctuation (AA 2 
as a function of pressure p for different ring sizes N. 



6.1 Numerical Scheme 




Fig. 4. Schematic view of a section of the simulated polygon. 
A randomly chosen vertex (marked by a circle) can be moved 
to a single position only (marked by a diamond) so as to main- 
tain the lengths of the two links attached to it constant. The 
area enclosed by the resulting rhomb is the difference in total 
polygon area for the given step. 



We consider a polygon of N equal edges, defined by the 
2D coordinates of its vertices. An off-lattice simulation 
is used, i.e., the positions of the vertices are defined in 
continuous space. At each step a random vertex is chosen 
and moved to the only other position that satisfies the 
edge- length constraint (Fig. The difference in energy 
between the two steps is proportional to the difference in 
total area, which in turn is simply given by the area of 
the rhomb composed of the two edges prior to and after 
the move (see Fig.0}. This way each step takes only 0(1) 
operations. The move is subsequently accepted or rejected 
according to the Metropolis criterion. 

The initial configuration is a stretched, regular poly- 
gon. This initial condition and the dynamics defined above 
imply that the polygon angles are restricted to change in 
discrete quanta of ±2w/N. Thus, although the algorithm 
is off-lattice, the simulated ring is a discrete variant of a 
freely jointed chain which strictly coincides with the freely 
jointed model only for N — > oo. 

The simulations were performed for N between 50 and 
3200, and for p between and 4p c . Away from the crit- 
ical point, the number of steps required for equilibration 
is 0(N 3 ), but near p c the simulation length must be ex- 
tended due to critical slowing down |21| . This limited our 
investigation of the transition to N < 3000. 



6 Monte Carlo Simulations 

Since our Flory argument and mean-field theory may fail 
near the critical point, we conducted Monte Carlo simu- 
lations to obtain the mean area (A), mean-square radius 



6.2 Results 

Figure shows simulation results for the mean area as 
a function of pressure for different values of N. When p 
is scaled by p c ~ A^" 1 and (A) by A^ (below p c ) or by 
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Fig. 5. Mean area as a function of pressure below (A) and 
above (B) the critical point. The pressure is scaled by p c = 



4n/N, and the area by N in (A) and by A n 



JV7(47r) 



(B). Symbols show the results of MC simulations for different 
values of N. The dashed lines show the prediction of the mean- 
field theory in the limit N — > oo. 

^4max ~ N 2 (above p c ), the data below and above the 
transition collapse onto two universal curves, thus con- 
firming the predicted scaling laws for the crumpled and 
smooth states (see Sections and 0J. However, while the 
data well below p c coincide with the scaling function ob- 
tained from the mean-field approximation, the data above 
the critical pressure collapse onto a different curve. 

The simulation results for the mean area and com- 
pressibility at p = p c as a function of N are shown in Fig. 
EK- The reduced compressibility, defined in Eq. jltty. was 
calculated from the measured mean-square area fluctua- 
tion as k — (4:Tt/N)(AA 2 )/(A). We get 



(A(p = Pc )) = (0.102 ± 0.007)iV 149±001 
«(P = Pc) = (0.56 ± 0.09)Ar°- 495±0 - 025 . 



(28) 
(29) 



Hence, the mean-field exponents for the scaling with N, 
Eqs. I|16(l and (|19|> . are confirmed. The predicted prefactor 



of the compressibility is within the standard error of the 
fit while that of the mean area differs by 4 standard errors. 

We also measured from simulations the dependence of 
the mean-square radius of gyration on N, as depicted in 
Fig. As predicted, below the critical pressure we find 
(R 2 g ) ~ N uol±om , and above it ^ _ jvi.985io.on_ At 

P = Pc we get 



(R 2 (p = p c )) = (0.043 ± 0.025)7V 



1.46±0.13 



(30) 
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Fig. 6. (A) Mean area and compressibility at p = p c as a 
function of N obtained by MC simulations. The fits, (A(p c )) = 
(0.102±0.007)A 149±0 - 01 and K (p c ) = (0.56±0.09)A°- 495±0 '° 25 , 
are given by the solid lines. (B) Mean-square radius of gy- 
ration as a function of N at different values of p — \p c , Pc 
and 2p c with the best fits (solid lines) (i? 2 (±p c )> = (0.12 ± 
0.01)Ar 101±0 ' 02 , (RliPc)} = (0.043 ± 0.025)A' 1 - 46±0 ' 13 , and 
{R 2 g {2p c )} = (0.019 ± 0.001) A- 1!)85±0 011 . 



For illustration we show in Fig.|7|four randomly chosen 
conformations of an 1600-segment ring at the critical pres- 
sure. The shapes vary in size significantly due to critical 
fluctuations. 



Fig. 7. Four random conformations of a ring with N = 1600 at the critical state as obtained by MC simulations, demonstrating 
the critical fluctuations. (The positions of the rings have no significance.) The spacing between gridlines is ten times the link 
length. The dotted circle, shown for reference, has an area equal to the mean area of the ring at this state. 



The simulation results for the smooth state, p > p c , 
are shown in Fig. |SJ where they can be compared with the 
transfer-matrix calculation, Eq. (|26[) . and the mean- field 
result, Eqs. (|10|l and Qllfl. On the one hand, there is good 
agreement with the transfer-matrix calculation for p > 5, 
particularly compared to the mean- field result. On the 
other hand, the mean-field theory succeeds in reproducing 
the phase transition, whereas the zeroth-order transfer- 
matrix calculation is invalid for these low pressures. 




5 10 15 20 

A 
P 



Fig. 8. Mean area (scaled by A max ) as a function of pressure 
(scaled by p c ), as obtained by the zeroth-order transfer-matrix 
calculation (solid curve), mean-field approximation (dashed 
curve), and MC simulations for N = 1600 (error bars). 



6.3 Critical Exponents 

We have confirmed the predicted exponents relating (A), 
(Rg) and k at p = p c with the system size N. We now turn 
to the exponents characterizing the divergence with \p— 1|. 
As in any continuous phase transition, the critical fluctu- 
ations make it hard to accurately measure these critical 
exponents from simulations. Instead, we choose the more 
reliable route of finite-size scaling to obtain the relations 
between them. 

Let us divide the ring at the critical state into Pincus 
correlation blobs |17II22| of g monomers and diameter £ 
each, such that within a blob the polymer behaves as an 
unperturbed random chain, i.e., g ~ £ 2 , whereas the chain 
of blobs is stretched by the pressure. The perimeter of the 
ring is R ~ (N/g)£ ~ N/£. We have already established 
that the perimeter of the ring at p — p c scales as R ~ 
iV 3 / 4 . (See Fig. 03.) Thus, for a finite-size system, we get 
a correlation length which scales at p — p c as 

Z(p=p c )~N 1 / 4 . (31) 

On the other hand, close to the critical point the correla- 
tion length diverges as £ ~ \p — 1|~", the compressibility 
as k ~ \p — 1| -7 , and the order parameter increases as 
M ~ \p — Using Eq. I|31|l and the numerically estab- 
lished results, M = R/N - iV- 1 /4 an( j K „ N i/2^ we 
obtain the relations 

P = v, 1 = 2v, (32) 

which hold for the mean-field universality class. (Note that 
the exponents a and S are irrelevant for this system, since 
both the ordering field and temperature are incorporated 
in the single parameter p.) 
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7 Discussion 

We have demonstrated that the swelling of a 2D freely 
jointed ring due to a pressure differential exhibits a second- 
order smoothening transition of the mean-field univer- 
sality class. Below the critical pressure the ring behaves 
as a random walk, with both the mean area (A) and 
mean-square radius of gyration (-Rg) proportional to N. 
In this crumpled state the mean area obeys the scaling 
(A) = Nff(p/p c ). Mean-field theory accurately captures 
the scaling law as well as the scaling function f^. (See 
Fig. EJA.) This lies in the fact that, for an unstretchcd 
chain, the Gaussian-spring description in the mean-field 
calculation and the actual freely jointed model coincide 
as N — ► oo. 

Above the critical pressure (A) and (Rg) are propor- 
tional to TV 2 . For this smooth state we have found a new 
scaling behavior, (A) = N 2 f^ip/Pc) (Fig.[5j3). The mean- 
field theory correctly gives the scaling law, yet fails to pre- 
dict the correct scaling function f^. This is because in a 
stretched state the entropy of a chain of Gaussian, vari- 
able springs is much larger than that of a freely jointed 
chain of rigid links. For p ^> p c we have calculated 
exactly [Eq. and Fig. EI - 

At the transition point itself a new state of polymer 
chains has been discovered, with (A) and (R 2 ) propor- 
tional to N 2 " c , v c = 3/4. The swelling exponent v c turns 
out to be identical to that of a 2D self-avoiding walk, al- 
though the physical origins of the two exponents are un- 
related; the freely jointed ring in this intermediate state 
between crumpled and smooth behaviors contains numer- 
ous intersections, as is clearly seen in Fig. {7\ 

Our discussion has been restricted to freely jointed, 
self-intersecting chains. We have commented in Section 
that the addition of self-avoidance to the Flory free en- 
ergy removes the second-order transition. Thus, the cur- 
rent work suggests a different perspective for the previ- 
ously studied swelling behavior of 2D self-avoiding rings 
6 , 7 . The large-inflation scaling regime thoroughly dis- 
cussed in those works can be viewed as a broadening of 
a phase transition present in the non-self-avoiding model, 
the self-avoidance thus acting as a relevant parameter for 
the critical point discussed in the current work. 

One of the strongest assumptions underlying our anal- 
ysis, as well as those of Refs. |51ll0llll| . is the replacement 
of the actual geometrical area of the ring with its alge- 
braic area. An important issue is how this assumption 
affects our results concerning the transition. It is clear 
that negative contributions to the algebraic area are sig- 
nificant in the crumpled, random-walk state and insignif- 
icant in the smooth state. The key question, therefore, is 
whether they are statistically significant at the transition. 
Returning to Fig. and the blob analysis presented in 
Section 16.31 we infer that the negative area contributions 
lie within the correlation blobs. We have shown that each 
blob contains g ~ £ 2 ~ N 1 / 2 monomers. The deviation of 
the geometrical area of a blob from its algebraic area is 
(a 2 ) 1/2 ~ £ 2 ~ N 1 / 2 . (Recall that the blob contains an 
unperturbed chain with zero mean algebraic area.) The 



ring contains n = N/g ~ jV 1 / 2 such blobs. Hence, the 
total deviation of the geometrical area from the algebraic 
one is (n(a 2 }) 1/2 ~ N 3 / 4 . In the limit N — >• oo this is 
negligible compared to the mean area of the ring at the 
critical state, (A) ~ iV 3 / 2 . Thus, we conjecture that the 
same smoothening transition as the one reported here will 
be found also in a model which considers the geometrical 
area rather than the algebraic one. 

We thank M. Kroyter, H. Orland, P. Pincus and T. Witten 
for helpful discussions. This work was supported by the US- 
Israel Binational Science Foundation (Grant no. 2002-271). 



Appendix: 2D Freely Jointed Chain under Ten- 
sion 

In this Appendix we recall the results for the partition function 
of a 2D, open, freely jointed chain subject to a tensile force. 
These results are used in Sections and 

Consider a 2D freely jointed chain composed of N links of 
length I = 1. A chain configuration is defined by N 2D unit 
vectors, {u., }■,=!,. .^jv, specifying the orientations of the links. 
One end of the chain is held fixed while the other is pulled by 
a force f (measured in units of ksT/l). The partition function 
of the chain is 



Z(f,N)= / ijdaje'^ 



3=1 

2n X W 



(A.l) 



/ r2ir \ JV 

= [J d8e fcos0 ) = [2nI (f)] N . 
The mean end-to-end vector is obtained from Eq. 1A.11 as 



R= V f \nZ = N^Q-f. 

M/J 



(A.2) 



The mean end-to-end distance in the limit of weak force, to 
two leading orders, is 



R/N ~ f/2 - /716, 
leading, upon inversion, to 

/ ~ 2R/N + (R/N) 3 . 



(A.3) 



(A.4) 



Finally, the free energy for fixed end-to-end distance (to two 
leading orders in small R) is 

F(fl) = -lnZ[f(R)]+f.R^^ + ^. (A.5) 

This yields the usual Gaussian term, F e i, and the first correc- 
tion due to inextensibility, -Finext, used in Eq. @. 
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